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Quantum dots are nanoscopic systems, where carriers are confined in all three spatial directions. 
Such nanoscopic systems are suitable for fundamental studies of quantum mechanics and are can¬ 
didates for applications such as quantum information processing. It was also proposed that linear 
arrangements of quantum dots could be used as quantum cascade laser. In this work we study the 
impact of electron-electron interactions on transport in a spinful serial triple quantum dot system 
weakly coupled to two leads. We find that due to electron-electron scattering processes the transport 
is enabled beyond the common single-particle transmission channels. This shows that the scenario 
in the serial quantum dots intrinsically deviates from layered structures such as quantum cascade 
lasers, where the presence of well-defined single-particle resonances between neighboring levels are 
crucial for device operation. Additionally, we check the validity of the Pauli master equation by 
comparing it with the first-order von Neumann approach. Here we demonstrate that coherences are 
of relevance if the energy spacing of the eigenstates is smaller than the lead transition rate multiplied 
by h. 


I. INTRODUCTION 

Electron-electron interaction effects in quantum dots 
have been a topic of active research within the past 
decades.The ability to confine a finite number of 
charged particles with advanced fabrication techniques in 
these structures opened up a possibility for testing dif¬ 
ferent physical theories such as charge and conductance 
quantization,Coulomb blockade,exciton forma¬ 
tion, mention a few. 

Here we focus on the electric transport through a se¬ 
rial arrangement of multiple quantum dots. Experimen¬ 
tally they can be realized in different ways and prominent 
examples are the gating of a two-dimensional electron 
gas,^^ cleaved edge overgrowth structures,stacked 
self-organized quantum dots,^^ nanowires either with ex¬ 
ternal gates^^’^^ or embedded heterostructures,^^ or the 
arrangement of atoms by a scanning tunneling micro¬ 
scope.^^ The restriction of phase space in such low¬ 
dimensional structures is reducing the scattering rates 
substantially and therefore such structures have been 
suggested for a wide range of applications ranging from 
quantum information processing^^’^^ to quantum cascade 
lasers. 

Electron transport through these structures has been 
widely used for level spectroscopy.Generally, one as¬ 
sumes, that the transport through quantum dot systems 
is dominated by specific resonances. These occur due 
to the alignment of energy levels in individual dots with 
those in neighboring dots as well as with the chemical 
potentials of metallic leads. This provides specific con¬ 
ditions for transport, which are resolved as current or 
conductance peaks for varying external parameters, such 
as the voltages at different gates. Such resonances may 
even refer to states with different energies due to the 
emission of phonons with a specified frequency^^’^^ or 
Auger processes.But even in this case the existence of 
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Fig. 1. Sketch of the serial triple dot structure attached to the 
source and drain leads. The electrons from the source/drain 
lead can tunnel into the left/right dot with the rate Tl/^r- 
The levels within the dots are coupled by the tunneling ma¬ 
trix Q. Inset: Current as a function of the left Fermi level 
taking into account only the Coulomb scattering term Use- 
Parameters as in Table I. 


specific resonances is the guiding theme of studying mul¬ 
tiple dot systems. However, with increasing number of 
dots the number of resonance conditions becomes large 
and difficult to satisfy simultaneously. Taking into ac¬ 
count growth imperfections as well as undefined locations 
of impurities with fluctuating charges, a strong suppres¬ 
sion of current is expected with an increasing number of 
dots.^^’^^ 

Electron-electron interaction is naturally occurring in 
all electronic devices and affects transport both by scat¬ 
tering (such as the Auger term) and level shifts. Eor 
systems with many degrees of freedom, such as bulk or 
layered systems, the continuum of states justifies usually 
a mean-field description, so that one can apply effective 
single-particle levels with renormalized energies. In this 
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case the resonances occur for different parameters, but 
the essential principle remains. A very successful exam¬ 
ple for this concept are quantum cascade lasers, whose 
operation is based on a clever design for the alignment of 
such single-particle levels.For quantum-dot systems, 
however, any mean-field model is very questionable, as 
one replaces the interaction between quantized charges 
by the interaction of a charge with an averaged quantity. 

This work analyzes finite bias electronic transport in 
both spinless and spinful triple quantum dots coupled in 
a serial configuration. The considered system serves as a 
simple model for longer arrangements such as quantum 
dot superlattices^^ or possible dot-based quantum cas¬ 
cade lasers (QCLs).^^’^^ The main interest of the present 
work is to address the influence of different parts of the 
electron-electron (ee) interaction inside the triple dot on 
the electrical current. Triple quantum dots in other kind 
of arrangements like triangular shape was extensively 
studied both theoretically and experimentally.^^’^^’^^”^^ 
We show that the Coulomb interaction between electrons 
opens up a large variety of different channels, which go 
far beyond the simple pictures of a few resonances, espe¬ 
cially when the spin degeneracy of the levels is included. 
We identify two main causes: (i) Coulomb scattering pro¬ 
vides a generic possibility for energy relaxation within 
the dot system, (ii) The multitude of many- particle 
states provides an enormous amount of different possible 
current paths. Additionally, this multitude of possible 
current paths give the possibility to study the applicabil¬ 
ity of simple Pauli master equation approach in realistic 
physical setup. In this work, we contrast the transport 
results obtained by Pauli master equation and first-order 
von Neumann approach. 


II. METHODS 

The system under consideration consists of three se¬ 
rial quantum dots sandwiched between metallic leads (see 
Fig. 1). In this section we specify the Hamiltonian, the 
different approaches for transport, and the parameters 
used in our calculations. 


A. Hamiltonian 


L^R stands for the left or right lead, a =t,i denotes 
the spin of the electron, and Ek£ = Ek E is the single 
particle energy of the electron in state k. The dispersion 
Ek in the lead is shifted by the chemical potential jui of 
the respective lead i. Also it is assumed that the lead 
states constitute a continuum with Ek G [—D^D] having 
a bandwidth 2D and a constant density of states up at 
the Fermi level. The dots are coupled to the leads by the 
tunneling Hamiltonian: 

+ h.c.) , (3) 

n,k£a 


with being electron creation operator in the dots, 
where n G {1,2, 3,4, 5} labels single-particle dot states 
as depicted in Fig. 1. Only the levels of the left {n = 
1,2) and right {n = 5) dots are directly coupled to the 
left and right lead, respectively. Additionally, we assume 
that all couplings to the leads have the same magnitude 
and phase. This means that we have the following non¬ 
vanishing tunneling amplitudes 


tlL=t2L=tL, t5R = tR, (4) 

which are parameterized by the tunneling rates ^^l/r = 
^^F\tL/R\‘^‘ We note that for symmetric structures a 
sixth level in the third dot below E^ is expected. How¬ 
ever, it is not considered here, as it hardly contributes to 
the current flow. 

The Hamiltonian of the dots is given by 


-Hd — ^ ^ ^nd\i(jdn(j ^ ^ ^ ^nmd\i(j^m(j “1“ -^ee? (5) 

ncr nm<j 


where ^nm describes the coupling between states n and 
m in neighboring dots. The electron-electron (ee) inter¬ 
action is described by the Coulomb Hamiltonian 




1 

2 


mnkl 

/ 


(6) 


For a system that has more than one confined electron, 
this part plays an important role and we will discuss the 
different matrix elements Vmnki Ihe following subsec¬ 
tion. In general the Coulomb matrix elements read 


The serial triple dot is modelled by the following 
Hamiltonian: 


= + + (1) 

where 

^LR = ^ ^ki^kia^kia^ (^) 

kia 

describes the source and drain leads as reservoirs with a 
continuum of noninteracting electrons, where denote 

the electron creation operators in the leads. Here £ = 


^mnkl — 


ATTSrSo 


d'^r 




r — r 


( 7 ) 

where is the spatial part of the single particle state 

m, Sr and Sq are the relative and vacuum permittivity, re¬ 
spectively. We neglect all terms, where either m and I or 
k and n belong to different quantum dots, as their over¬ 
lap would be vanishingly small. Furthermore, terms con¬ 
necting levels of next-nearest neighboring dots are small 
and neglected as well. The remaining terms can be cat¬ 
egorized into Intradot and Interdot interactions and are 
separately treated below. 
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1 . Intradot Interaction 


For intradot interaction all the levels mnkl are consid¬ 
ered to be in the same dot. By employing the normaliza¬ 
tion condition for the wave function, the direct elements 
can be estimated as: 




AnerSocr 


= U, 


(8) 


where a = ((r — (r))^) is the standard deviation for the 

spatial extension of the dot wave functions. Another set 
of interaction matrix elements that has to be taken into 
account are Vmnmn (for n 7 ^ m), which act as exchange 
terms for equal spins and scattering terms for different 
spins. Trying different test wave functions, we observe 

^mnmn ~ Uqx with Uex — • (9) 


2 . Interdot Interaction 

The direct interaction between two states in the neigh¬ 
boring dot can be approximated in a similar way as 
Eq. ( 8 ): 


Vm 


AirerSod 


= Un, 


( 10 ) 


where d is the distance between the centers of the dots. 
The terms with different combinations of indices, are es¬ 
timated by a Taylor expansion of l/|r — r'| around the 
centers of the respective dots R^, Rj, see Ref. [45]. Using 
|Ri — Rj| = d for neighboring dots, we find 


o 

II 

u = 10 

Fz. = 0.1 

Qi3 = -0.05 

o 

II 

(M 

Uex = 2 

Fr = 0.1 

Qi4 = 0.1 

Es = 20 

Un=S 

fiL = 50 

U 23 = 0.1 

o 

II 

Use = -0.2 

/an = 10 

U 24 = 0.2 

E5 = 20 

0 

1 

II 

0 

A;bT = 1 
D = 10^ 

U 35 = 0.1 

U 45 = 0.2 


Table I. Parameters used in the calculations if not mentioned 
otherwise. All energies are in meV. 


40 nm) embedded between InP barriers (thickness 3 nm). 
Similar structures have been recently fabricated. 

The values of the energies Ei and couplings ftnm are 
estimated by a tight-binding superlattice model (see 
Ref. [48]), as outlined in the Appendix C. This model also 
provides the dipole matrix element |s 2 i| =8 nm, which 
points in the direction of the nanowire, and is used to es¬ 
timate the Coulomb matrix elements as outlined above. 
Furthermore we use a = 11 nm and d = 43 nm as well 
as Sr^inAs = 15.^^ In addition, the Coulomb matrix el¬ 
ements were calculated numerically from Eq. (7) using 
wave functions in cuboids (40 nm x 35 nm x 33 nm) 
representing each dot. This result was in good agree¬ 
ment (about 10 % deviations) with the approximations 
addressed above. As we are only interested in general 
features, rather than in modelling a specific device, we 
used rounded values for all quantities in order to allow 
for an easy recognition of scales in the plots. The specific 
values are given in Table I. 


C. Transport Calculation 


( 11 ) 

AirSrSo d'^ 

Vmnkl ~ ^ ^ 

with the intradot dipole matrix element 

Snm = J d^rip’^{T^)ripm{r). (13) 

The Udc and Use terms can be interpreted as a dipole- 
charge interaction and dipole-dipole scattering terms, re¬ 
spectively. The Use term is responsible for the Auger 
process and is crucial for the current flow in our system. 
The sign of Ude fo Eq. (11), depends on whether the 
charge is on the right or left side of the dipole. 


B. Parameter values 

In order to obtain realistic and consistent parameters, 
we consider a particular nanostructure which is made 
from a nanowire containing three InAs wells (thickness 


We are interested in obtaining the current from the 
left to the right lead for considerably large bias. In order 
to do that, we diagonalize the Hamiltonian (5) and 
get the many-particle eigenstates |a), | 6 ),... of the triple 
dot. Expressed in this many-particle basis the tunneling 
Hamiltonian Ht (3) becomes 

Ht= Yi {TUi<y)\b){a\c,,^ + h.c.), (14) 

ah,k£a 

Tbai^cr) = ( 15 ) 

n 

Here we used the letter eonvention: if more than one 
state enters an equation, then the position of the let¬ 
ter in the alphabet follows the particle number (for ex¬ 
ample Ni, = Na E 1, Nc = Na E 2 , Na' = Na). In 
such a way the sum restricts to those combinations, 
where Nc = 1. To obtain the current through the 

device we use the first-order von Neumann (IvN) ap¬ 
proach,where coherent effects are included, and a sim¬ 
pler Pauli master equation where only populations are 
considered.The derivation of the governing equa¬ 
tions is presented in Appendix A and Appendix B. The 
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IvN approach treats the reduced density matrix in low¬ 
est order with respect to It is conceptually simi¬ 

lar to the Wangsness-Bloch-Redfield^^ (also presented in 
Appendix A), albeit, the Markov limit is done in a dif¬ 
ferent way. For our calculations we did not observe any 
differences in the current between these approaches as 
long as principal part integrals are neglected. We choose 
the temperature of the leads to be larger than the lead 
tunneling rate, ksT ^ F, which justifies the neglect 
of higher-order tunneling processes in the IvN and the 
Wangsness-Bloch-Redfield approach. Additionally, this 
suppresses any kind of Kondo correlations (dominating 
below the Kondo temperature Tk)^ as ksT^ < F.^^ 
Note that these approaches take into account electron- 
electron scattering by using the eigenstates of the dot- 
Hamiltonian. The Auger term of the Coulomb interac¬ 
tion with matrix element Use couples different configura¬ 
tions of many-particle states, which establishes connec¬ 
tions between the leads. This does not require any relax¬ 
ation terms inside the quantum-dot structure. Further 
relaxation processes, e. g., due to phonon-scattering, are 
entirely neglected here. If the emission of optical phonons 
is energetically not allowed, phonon scattering rates be¬ 
tween quantum dot states are of the order of 1 /ns or even 
smaller.Such a scattering process can provide back¬ 
ground currents of at most ^0.1 nA. As this is smaller 
than the current peaks addressed in the subsequent sec¬ 
tion, it is justified to neglect these processes in our study. 


III. RESULTS 

In the following we show calculations for varying levels 
Es and F^ 4 . This is motivated as follows: even in a nom¬ 
inally identical triple-dot system, growth imperfections, 
charged impurities and the electrostatic environment will 
modify the energy levels in a way which is difficult to pre¬ 
dict. We fix the excitation energy in the left dot, E 2 —Ei^ 
as a reference point and the level E^ of the right dot can 
be shifted by a source drain bias. However, the remain¬ 
ing levels Es and E^ are less controllable, albeit there is 
a possibility of gating. Thus it is a central question of 
practical interest, for which range of parameters 
current flow occurs. 

For an applied bias, standard considerations of coher¬ 
ent transport, such as the transmission formalism, pre¬ 
dict current flow, if there is a state connecting both leads. 
If we neglect the Coulomb interaction at all, there is only 
a very small current unless three levels align in a row 
(e.g. for El ^ E 4 ^ E^, 8 i case we do not consider). For 
the situation in Fig. 1 one would thus not expect any elas¬ 
tic transport through the structure. Neglecting phonon 
scattering, the current would thus be zero, even if the 
energies E^ and E 4 are varied. However, in Ref. [57] it 
was shown, that for this situation, electron-electron scat¬ 
tering enables a current flow: if the levels 1 and 4 are 
occupied, an energy-conserving scattering event creates 
the simultaneous transitions 1^2 and 4^3 (for E^-E^ 


= E 2 -E 1 ), with subsequent tunneling to the right lead 
via the state 5. The inset of Fig. 1 shows the current 
as a function of left Fermi level. As it is shown the cur¬ 
rent is decreasing drastically when the left Fermi level 
comes close to the excited state of the first dot. Then the 
level E 2 becomes occupied, and consequently the Auger 
scattering process is prevented due to the Pauli blocking. 
This is a distinct feature of the Auger driven transport in 
multiple dots, which discriminates it from other possible 
scattering processes. Thus, scattering-assisted processes 
dominate. 

A finite Coulomb scattering matrix element V 2341 = 
Use allows for an Auger process, where one electron re¬ 
laxes from level 4 to level 3, while transferring its energy 
to an electron being excited from level 1 to level 2. There¬ 
fore Use = — 0.2 meV is used in all calculations, while the 
impact of the other Coulomb terms is neglected in some 
calculations in order to demonstrate their relevance. The 
left and right Fermi levels are fixed to be /ix = 50 meV, 
IJ.R = 10 meV so that the leads fill level 1 and empty the 
levels 2 and 5. The calculations using the Pauli master 
equation and the IvN approach are presented in Figs. 2 
and 3, respectively. 


A. Spin-polarized levels 

At first we restrict to spin-polarized levels in the quan¬ 
tum dot. This could be achieved by having a large mag¬ 
netic field in the system. Figures 2a and 3a show the 
current as a function of Es and F^ 4 , when just scatter¬ 
ing elements Use are included. Only for specific values 
of Es and E 4 the electrons are able to pass through the 
system. This reflects the conservation of energy both for 
the electron-electron scattering and for tunneling, i. e., 
Es = E^ = 20 meV and E 4 = Ei = 40 meV. Otherwise, 
the electron transport is blocked (or limited to a back¬ 
ground current below 0.1 nA due to phonon-scattering 
in real structures as discussed above). Thus, any slight 
changes in the geometry of the dots and the configura¬ 
tion of the energy levels would prevent the current flow. 
Such a selective situation would be optimal for devices, 
which rely on well defined transitions, such as quantum 
cascade lasers. 

When ee-interaction with all the matrix elements is 
considered, we see that the current can flow for a wider 
range of parameters, as shown in Figs. 2b and 3b. This 
scenario is reflected by the addition energies, where a par¬ 
ticle can tunnel into or out of the triple dot. In Fig. 4 we 
display the differences in energy between all possible two 
and three-particle states {AE axis). In order to specify 
their relevance for single-particle transitions, we plot the 
respective transition probabilities for electrons to enter 
from either lead, while changing the state of the system 
between these states. For the case restricting to Coulomb 
scattering, there are only three distinct energies, where 
lead electrons can enter, which correspond to the ener¬ 
gies of the levels 1,2, and 5. In contrast, a larger variety 
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Fig. 2. Current through the triple dot system calculated using the Pauli master equation approach, (a) Spin-polarized system 
with just scattering Use included [Eq. (12)]. (b) Spin-polarized system with all Coulomb matrix elements included, (c) Spinful 
system with just scattering Use- (d) Spinful system with all Coulomb matrix elements. Parameters as in Table I. 


of excitation energies is relevant if the full Coulomb in¬ 
teraction is taken into account. This is fully consistent 
with the differences between the current plots of Figs. 3a 
and Fig. 3b. 


B. Spinful levels 

In the absence of an external magnetic field each single¬ 
particle energy level is doubly degenerate due to the 
Kramer’s theorem. Let us first consider the case with 
just scattering t/gc elements present (Fig. 3c). In a sim¬ 
ple picture, one could expect, that the current doubles 
compared to the case where the levels are spin-polarized 
(Fig. 3a). The maximum of the current in Fig. 3c shows 
indeed an increase by a factor of two. Furthermore, 
Fig. 3c shows that the resonance faintly extends along 


the intersection of E^ = 20 and E/^ = 40 lines. However, 
in a spinless case this resonance is extending just along 
the = 40 line. 

In order to understand the reason of the different be¬ 
havior in Fig. 3, we study a line for fixed = 38 meV in 
Fig. 5 and corresponding eigenstates as indicated by the 
squares. For parameter values at the green (right) square 
significant current is observed for the spinful levels, but 
current is blocked for the spin-polarized levels. At the 
red (left) square current is blocked for both cases. 

The five one-particle eigenstates of the Hamiltonian 
responsible for transport are illustrated in the green di¬ 
agram in the lower diagrams of Fig. 5. There is a strong 
coupling between Ei and F ^4 as well as between E^ and 
E^ thus the superposition of each of these two states cre¬ 
ates two new states with two new energies, which are 
referred to as F^i, 4 , L^ 3 , 5 , and E’ 5 , 3 . For the pa- 
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Fig. 3. Same as Fig. 2 but with the current calculated using the IvN approach, (a) The current has a single distinct peak, 
where E 3 ^ E 5 and E 4 ~ Ei and the resonance extends along E 4 = 40 meV line, (b) Inclusion of all possible Coulomb 
elements opens up additional channels for transport leading to more current peaks, (c) The maximum current becomes larger 
by a factor of two compared to the spin-polarized case, however, now the resonance extends faintly along the E 4 = 40 meV and 
E 3 — 20 meV lines, (d) For spinful system additional Coulomb elements open very many channels, which leads to significant 
background current. Parameters as in Table I. 


rameters at the right (green) square the conservation of 
energy is satisfied, {E 2 — Ei^ 4 ) ^ (^ 1,4 — ^ 3 , 5 ): which 
allows for Coulomb scattering. However, the state Ei ^4 
appears as the initial state in both parts for an Auger 
process. Thus, the Pauli principle can only be satisfied 
if the level is spin degenerate. Due to this reason the 
electrons are able to transfer through the triple dot in 
a spinful system, but the current is blocked in a spin- 
polarized system. Spin is definitely more than a factor of 
two here. We note that this transport channel is rather 
sensitive to actual couplings to the leads. The current is 
reduced if the coupling to the leads is increased as can 
be seen from dashed (green) curve. The coherences are 
responsible for this surprising reduction as discussed in 
the next section. 


Including all ee-interaction terms for spinful system 
provides a wide variety of single particle excitations and 
consequently resonance conditions are easier to satisfy 
than for the spin-polarized case addressed above. Fig¬ 
ures 2d and 3d show the evaluated current and we observe 
a multitude of peaks as well as significant background 
current of ^ 0.15 nA for a large number of energy level 
combinations. As each peak relates to a different current 
path through the structure, it is very difficult to address 
or identify a specific transport path in an experiment. 
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Fig. 4. Addition energies (horizontal) for the two-particle 
states in the spin-polarized system. The vertical axis shows 
the respective coupling strength for electronic transitions from 
either lead. In the left panels a point is drawn for each transi¬ 
tion. In order to resolve multiple transitions, the right panel 
sums Lorentzians with full width at half maximum of 5 meV 
with a peak given by the points in the left panel. 




Fig. 5. Cuts of contour plots in Fig. 3c [solid (black)] and 
Fig. 3a [dotted (red)] showing increased current for spinful 
system due to Auger process. Dashed (green) and dashed 
dotted (blue) gives current when the coupling F is increased. 
The left and right square panels show energy spectrum for 
the one-particle states at two values of energy E^. The units 
of F in the legend are meV. 


C. Role of coherences: comparison of IvN and 
Pauli approaches 


In the following, we discuss the role of coherences, 
which in the IvN approach are defined as the off-diagonal 
elements of the reduced density matrix in the many-body 
eigenbasis of i^D (see Appendix A). The Pauli master 
equation approach neglects any such coherences, which 
are known to be of relevance, if /S.E < F, where /S.E is 
the separation between two relevant levels and F is the 
transition rate in units of energy. Typical examples for 
this situation have been already discussed in Refs. [58 
and 59]. In order to verify our results obtained by the 
Pauli master equation, we compare the IvN approach, 
which keeps the coherences of the system and takes into 
account the tunnel transitions to the leads in lowest or¬ 
der. If the temperature k^T of the leads surpasses the 
transition rate F, the IvN approach is believed to provide 
reliable results for the currents, as shown by comparison 
with higher-order approaches. 

By comparing Figs. 2 and 3 we see that if all interac¬ 
tions are included [(b) and (d)] the peak structure in the 
IvN and Pauli approaches is similar. However, the thin 
line of resonance for £^4 £^3 + 20 is substantially re¬ 

duced in the more advanced IvN simulation. Along this 
line the 2-particle state occupying the levels 1 and 3, and 
the 2-particle state occupying levels 4 and 5 are degener¬ 
ate. Via the second-order couplings they get mixed and 
anticross with a small splitting of A£ ^ 2 /icV. This en¬ 
ergy difference is much smaller than the transition rates 
F = 0.1, which explains, why this narrow line is mostly 
an artificial result in the Pauli master equation. The 
splitting of the 2 -particle states is enhanced to A£ ^0.1 
around £4 = 35, £3 = 15 for spinless and £4 = 32.5, 
£3 = 13 for spinful cases, where further 2-particle states 
are in resonance. This corresponds to a broader peak, 
which is also visible in the IvN simulation. If the cou¬ 
pling to the leads is increased, coherences become even 
more important, which leads to further reduction of these 
resonances. On the other hand, we note that the IvN cal¬ 
culations match the Pauli master equation result if the 
coupling to the leads F becomes vanishingly small com¬ 
pared to the energy A£ splittings between the many- 
particle states with the same number of particles, i. e. 
F/A£ 0. 

Similar considerations hold for the case when just scat¬ 
tering Use is included [see (a) and (c) in Figs. 2 and 3]: for 
spin-polarized levels (a) thin lines of resonances with high 
current appear at £4 ~ 40 and £4 « £3 + 20 (faintly) 
for the Pauli master equation, which gets suppressed by 
coherences. However, in the spinful case (c) the Pauli 
master equation provides a very different scenario com¬ 
pared to the IvN approach. In Fig. 2c many resonances 
are pronounced, which get completely diminished by the 
IvN approach in Fig. 3c, leaving just an extended res¬ 
onance along the intersection of £3 = 20 and £4 = 40 
lines. 


































D. Neglect of principal part terms 
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In the IvN approach calculations of Figs. 1-5 we have 
In the IvN approach calculations of Figs. 1-5 we have 
neglected the principal part integrals [see Eq. (A20)]. At 
the resonance position = 0 this integral gives a loga¬ 
rithmic divergence, which is cut-off by the temperature. 
So it could be expected that for large enough tempera¬ 
ture these terms should not be problematic. However, in 
Fig. 6 , we show a calculation for spin-polarized quantum 
dots, with all interactions and the principal part integrals 
included. We see that regions of negative current appear 
(flow against the bias), and at particular points the cur¬ 
rent also has divergencies (see Fig. 6 b). At these points 
of divergent current the positivity of the populations ^55 
is highly violated and the elements ^ 55 / acquire divergent 
structure as well. This unphysical behavior of the current 
motivates the neglect of the principal part integrals. 

Finally, the principal parts have the structure of the 
renormalization of the energy spectrum (Lamb shift^^) 
as can be seen from perturbation theory for energies^^ 
or real parts of self-energies in Green’s function formal¬ 
ism.^^ This suggest that this type of terms should be 
resummed to get an effective Hamiltonian in which new 
eigenspectrum of the quantum dots is obtained and the 
energy differences E 5 — Ea are renormalized.^^’^^ 


IV. DISCUSSION 
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The main results are collected in Fig. 3, which dis¬ 
plays the current through both spin-polarized and spin¬ 
ful triple quantum dot system, calculated using the IvN 
approach. For the single spin case restricting just to scat¬ 
tering Use, a single peak is seen in Fig. 3a, which is easily 
predicted by standard single particle states. Such reso¬ 
nances are well known and are commonly used for device 
design. However, taking into account spin degeneracy of 
the states and all interactions. Fig. 3d displays an en¬ 
tirely different picture, where current flow is spread over 
a wide range of parameters where unexpected paths be¬ 
come of relevance. As a consequence, designing devices 
based on multiple quantum dots is questionable, if one 
restricts to single-particle models even if the mean-field 
is included. For example, for dot-based quantum cascade 
lasers, it is difficult to accomplish the desired specificity 
of electron injection into the upper laser level. On the 
other hand, the multitude of channels creates significant 
background current for wide range of energy level con¬ 
figurations, which makes transport in multiple quantum- 
dot devices less sensitive to size fluctuations and random 
charges. 

We have compared the Pauli master equation ap¬ 
proach, which neglects coherences, and the IvN approach 
where all density matrix elements of the reduced quan¬ 
tum dot system are taken into consideration. It was 
shown that the qualitative picture of resonances is cor¬ 
rectly captured by the simpler Pauli master equation. 


Fig. 6. Current for spin-polarized system with all interactions 
calculated using the IvN approach [similar to Fig. 3b], when 
the principal part V integrals are included [see Eq. (A20)]. 
(a) Due to divergent current the values below —0.1 nA and 
above -1-0.3 nA are filtered, (b) Cut showing the comparison 
of calculated current with and without the principal parts. 


However, there are cases where the current is highly over¬ 
estimated by this approach and this can be seen by com¬ 
paring Eigs. 2c and 3c. These points are traced back to 
specific level configurations with the lead coupling larger 
than the energy splitting, i. e., F > AE, where coher¬ 
ences strongly reduce the current. This also allows for 
situations where the current drops with increasing lead 
coupling F as shown Fig. 5. Additionally, the relevance 
of principal part integrals in the IvN approach was ex¬ 
amined and it was shown that it can lead to highly un¬ 
physical results (current flowing against the bias), which 
suggests that neglecting such terms is reasonable. 
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Appendix A: First-order von Neumann (IvN) 
approach 

The dynamics of the full density matrix p of the system 
is described by the von Neumann equation: 

(Al) 

We solve the above equation approximately using the 
IvN approach.In this approach only the density 
matrix elements, which connect the states differing by 
just one electron or hole excitation, are considered. Ad¬ 
ditionally, such a treatment is valid if the temperature T 
is larger than the tunneling rates, T ^ T. In the other 
limit, T ^ r, different approaches based on numerical 
renormalization group,quantum Monte Carlo simula¬ 
tions,^^ or higher order expansions should be used.^^’^^“^^ 
The density matrix elements are defined as 

P^albg' = (“S'IpIV), (A2 ) 

where \bg) = \b) 0\g), with \b) denoting the eigenstate of 
the dot Hamiltonian (5) and \g) denoting the eigenstate 
of the lead Hamiltonian (2). Here the label n provides the 
number of electron or hole excitations needed to trans¬ 
form 1^) into 1^'). For example, we will consider matrix 
elements of the type 

Pblb'g = {b9\p\b'g), 

Pbl-K,ag = {^9 - I^\p\a9)- 
Here we have introduced the following notation 
K = k, cr; 

\bg + K) = | 6 ) ® 4 | 5 )> 

\bg-K) = \b)(S)cJg). 

By neglecting all the density matrix elements with more 
than one electron or hole excitation n > 2 from Eq. (Al) 
we obtain the equations 

+ TbaA^l)P^alg+.ub'g(-^)^^^ 

+T,,,(5i4cj5) 


(A3) 

(A4) 

(A5) 


Here denotes the number of electrons in the state \b). 
Note that all indices with subscript 1 like ai, ci, ki are 
summed over and the letter eonvention introduced in the 
main text is used. Also the matrix element (^|cjj.c^|^) 
in the second term of the right-hand side of Eq. (A7) 
corresponds to the requirement that there is an electron 
in the single particle state corresponding to k. Addi¬ 
tionally, phase factors like (—1)^^ appear due to order 
exchange of the lead operators with the dot operators, 
i.e., c^(\b) 0 14) = {-l)^^\b) <S)cJg). 

Summing Eqs. (A6) and (A7) over all the lead states 
\g) we get 

- (« i ) T ,, 6 '(« i ) - $ rA « i ) T „, 6 ,( Ki ) 


(A9) 


^ (Ec -E,- S6)$ib4 

where we introduced the following notation 

9 

/fc = fki = {e-xpiEk/k^Te] + 1)“\ 
f-n = 1 - fki- 


Here we have also assumed that the electrons in the leads 
are thermally distributed according to the Eermi-Dirac 
distribution / and that this distribution is not affected 
by the coupling to the quantum dots. This assumption 
leads to the following relations: 


(AlO) 


^Pbg,b'gi9\cicjg) sa 


EpI 


[ 0 ] 

^g—K,,b'g—K, 


f 


(All) 


Eor the stationary state we assume the conditions 

(A12) 

which allow to write in terms of as 




0 , ih|$W(4 = o, 


^ib4 = 


— Ec Ef, ip 


(A13) 


Here we have added a positive infinitesimal p = +0 to 
ensure a proper decay of initial conditions, which can be 
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seen by formally integrating Eq. (A9): 


X - $ E \( i ') rc ,&(«)/-.). 


(A14) 


After performing a Markov approximation in the above 
integral, ^ and setting t +oo we also 

obtain Eq. (A13). We note that there is another possi¬ 
bility to choose the time dependence of 

(A15) 

In this case we obtain, 


which appears after performing the k-sums using a flat 
density of states approximation, i. e., ^ up dE, 

with Up denoting the density of states at the Eermi level 
and 2D denoting the bandwidth of the leads. In our 
calculations we assume that the bandwidth of the leads is 
the largest energy scale. Note that, in the limit D —^ -Foo 
the results become bandwidth independent^^ and that is 
what we also have checked in our numerical simulations. 
We note that the IvN approach does not include any 
broadening effects of the quantum dot levels due to leads, 
which in principle could be included by different means 
(e.g., see Refs. [65, 70, 73, and 74]). 

Einally, we are interested in the current going from the 
lead £ into the quantum dots, which is given by 




-Eav ~ Ec + if] 

Ef^ — Ecl + Ef, + ir] 


(A16) 


r{t) = eJ 2 

kcr 

ka 


(A21) 


which together with Eq. (A8) provides the Wangsness- 
Bloch-Redfield approach.The denominators in 
Eq. (A16) differ from Eq. (A13) if non-diagonal density 
matrix elements are relevant. The oscillatory be¬ 
havior in Eq. (A15) is suggested by the first right-hand 
side term in Eq. (A8). 

After combining Eqs. (A8), (A12), and (A13) we get 
the IvN approach equations for the steady state: 


0 =^bb' {Eb — Ep) 

+ Y .^ bb "[ Y .^ i " a , ab'ka 

b"i a 

b"£ c 


E pi 

J- b"c,cb'-^cb 
c 

E pi j£-* 

J- ba,ab"-‘-b'a 
a 


aa'£ 


cc'£ 



(A17) 


Here A{t) = denotes the Heisenberg evolu¬ 

tion of an operator A. In the steady state the current is 
obtained from Eq. (A21) in terms of as 


I - 2 " 


im ^ 2^ i bcpb'-^-cb ^b'b ^ bcp'b^cb ^cd 


cb 


(A22) 

which is the main output presented in the paper. Eor 
all the calculations (except for Eig. 6) the principal part 
integral V in Eq. (A20) is neglected. The motivation for 
neglecting these terms is given in Section HID. 


Appendix B: Pauli master equation 

The Pauli master equation can be obtained from the 
IvN approach by neglecting the coherences ^bbd h ^ h'. 
In such a case for the populations Pb = ^bb we obtain 
the equations: 


Additionally, we impose the normalisation condition for 
diagonal density matrix elements: 

^^66 = 1- {E18) 

b 


E [PaAa^bfi+Pba) “ PbTLafi-Pba)] 

a£ 

+ Y.[PcTUd{-pi,) - PPUji+pi,)] = 0 , 


Here in Eq. (A17) the tunneling rate matrix T is defined 
as 

^ba,a'b' ~ E na{e<7)Ta'b'{e<7), (A19) 

a 

and the following integral was introduced 

2^4a = P r - inf{±pL)0{D - \pI\), 

J-D ^ - Pba 

Pba = Eb- Ea- Pi, (A20) 

f{E) = {eME/k^T] + 1)-!, 


where we have denoted d^b = dbM = = ^ba,ab- 

Using the populations Pb the steady states current is ex¬ 
pressed as 

T,[Pa^i^bf(+Pia) - PbrLafi-Pba)]- (B2) 

ab 

Appendix C: Motivation for parameter values 

We consider a superlattice structure with 40 nm InAs 
wells (meff = 0.026) and 3 nm InP barriers (meff = 0.08, 
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Fig. 7. Wannier levels for the superlattice structure at a volt¬ 
age drop of 21 meV per period. We use just three periods 
of this structure to model single particle levels in the serial 
triple dot. 


= 0-6 eV) providing the period d = 43 nm. 
The serial triple dot is described by just three peri¬ 
ods of such a structure. A standard calculation (see, 
for example, Ref. [75]) provides the minibands Ej^{q) « 


Ey + 2Ty cos{qd) with Ea = 7 meV, Ta = —0.058 meV, 
E}j = 28 meV, and = —0.23 meV. Figure 7 displays 
the corresponding Wannier functions, where an electric 
field of 21meV/(i is applied, so that the ground state a 
of one well is in resonance with the excited state b of the 
neighboring well. The tunnel couplings between Wannier 
states from the same band u in neighboring periods is Ty. 
Thus we identify 

• Qis = Ta = -0.058 meV 

• 1124 = II 45 =Ti) = 0.23 meV 

Furthermore, we obtain numerically the dipole matrix 
elements Snm = Smn = /between the 
states. The non-vanishing values are 521 = 543 = 8.16 
nm (which is also used for the calculation of Coulomb 
matrix elements for the interdot interaction after making 
a dipole expansion, see Ref. [57]) and S 41 = S 32 = S 53 = 
—0.21 nm. For the operation bias this provides the tunnel 
couplings between states from different bands Unm = 
—eESnm with the values 

• II 41 = II 53 = II 32 — 0.103 meV 

We note, that the same procedure is used for simulating 
quantum cascade lasers, where we obtain quantitative 
agreement for a variety of samples. 
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